Larger particle size distribution of environmental RNA compared to environmental DNA: a case study targeting the mitochondrial cytochrome b gene in zebrafish (Danio rerio) using experimental aquariums.

Environmental RNA (eRNA) analysis is conventionally expected to infer physiological information about organisms within their ecosystems, whereas environmental DNA (eDNA) analysis only infers their presence and abundance. Despite the promise of eRNA application, basic research on eRNA characteristics and dynamics is limited. The present study conducted aquarium experiments using zebrafish (Danio rerio) to estimate the particle size distribution (PSD) of eRNA in order to better understand the persistence state of eRNA particles. Rearing water samples were sequentially filtered using different pore-size filters, and the resulting size-fractioned mitochondrial cytochrome b (CytB) eDNA and eRNA data were modeled with the Weibull complementary cumulative distribution function (CCDF) to estimate the parameters characterizing the PSDs. It was revealed that the scale parameter (α) was significantly higher (i.e., the mean particle size was larger) for eRNA than eDNA, while the shape parameter (ß) was not significantly different between them. This result supports the hypothesis that most eRNA particles are likely in a protected, intra-cellular state, which mitigates eRNA degradation in water. Moreover, these findings also imply the heterogeneous dispersion of eRNA relative to eDNA and suggest an efficient method of eRNA collection using a larger pore-size filter. Further studies on the characteristics and dynamics of eRNA particles should be pursued in the future.

Comparative evaluation for the performance of environmental DNA and RNA analyses targeting mitochondrial and nuclear genes from ayu (Plecoglossus altivelis).

Environmental DNA and RNA (eDNA and eRNA; collectively eNA) analyses have the potential for non-invasive and cost-efficient biomonitoring compared with traditional capture-based surveys. Although various types of eNA particles, including not only mitochondrial eDNA but also nuclear eDNA and their transcripts, are present in the water, performances of eNA detection and quantification have not yet been evaluated sufficiently across multiple mitochondrial and nuclear genes. We conducted a tank experiment with ayu (Plecoglossus altivelis) to compare the detection sensitivity, yields per water sample, and quantification variability between replicates of each type of eNAs. The assay targeting the multi-copy nuclear gene exhibited a higher sensitivity than the assay targeting the mitochondrial gene, and both the target eDNA and eRNA concentrations per water sample were higher for the nuclear gene. On the contrary, variation in eRNA quantifications per sample does not necessarily correspond to that in eDNA, and the intra-sample quantification variability (represented as the CVs between PCR replicates) tended to be larger for eRNA than eDNA. Our results suggested that, even if suitable to the sensitive detection of species occurrence, the use of eRNA particularly derived from multi-copy nuclear gene may not be necessarily appropriate for the reliable assessment of species abundance. The findings in this study would help optimize eNA analyses for making biomonitoring and stock assessment in aquatic environments more efficient and reliable.

Validating post-enrichment steps in environmental RNA analysis for improving its availability from water samples.

Environmental RNA (eRNA) analysis is expected to inclusively provide the physiological information of a population and community without individual sampling, having the potential for the improved monitoring of biodiversity and ecosystem function. Protocol development for maximizing eRNA availability is crucial to interpret its detection and quantification results with high accuracy and reliability, but the methodological validation and improvement of eRNA collection and processing methods are scarce. In this study, the technical steps after eRNA extraction, including genomic DNA (gDNA) removal and reverse transcription, were focused on and their performances were compared by zebrafish (Danio rerio) aquarium experiments. Additionally, this study also focused on the eRNA quantification variabilities between replicates and compared them between the PCR and sample levels. Results showed that (i) there was a trade-off between gDNA removal approaches and eRNA yields and an excess gDNA removal could lead to the false-negative eRNA detection, (ii) the use of the gene-specific primers for reverse transcription could increase the eRNA yields for multiple mitochondrial and nuclear genes compared with the random hexamer primers, and (iii) the coefficient of variation (CV) values of eRNA quantifications between PCR replicates were substantially lower for those between samples. Including the study, further knowledge for the sensitive and precise detection of macro-organismal eRNA should be needed for increasing the reliability and robustness of eRNA-based biomonitoring.

Pooling of intra-site measurements inflates variability of the correlation between environmental DNA concentration and organism abundance.

Environmental DNA (eDNA) analysis can promote efficient ecosystem monitoring and resource management. However, limited knowledge of the factors affecting the relationship between eDNA concentration and organism abundance causes uncertainty in relative abundance estimates based on eDNA concentration. Pooling of data points obtained from multiple locations within a site has been used to mitigate intra-site variation in eDNA and abundance estimates, but decreases the sample size used for estimating the relationship. I here assessed how the pooling of intra-site measurements of eDNA concentration and organism abundance impacted the reliability of the correlative relationship between eDNA concentration and organism abundance. Mathematical models were developed to simulate measurements of eDNA concentrations and organism abundances from multiple locations in a given survey site, and the CVs (coefficient of variability) of the correlations were compared depending on whether data points from different locations were individually treated or pooled. Although the mean and median values of the correlation coefficients were similar between the scenarios, the CVs of the simulated correlations were substantially higher under the pooled scenario than the individual scenario. Additionally, I re-analyzed two empirical studies conducted in lakes, both showing higher CVs of the correlations by pooling intra-site measurements. This study suggests that it would make eDNA-based abundance estimation more reliable and reproducible to individually analyze target eDNA concentrations and organism abundance estimates.

Methodological considerations for aqueous environmental RNA collection, preservation, and extraction.

Environmental RNA (eRNA) analysis is expected to infer species' physiological information (health status, developmental stage, and environmental stress response) and their distribution and composition more correctly than environmental DNA (eDNA) analysis. With the prospect of such eRNA applications, there is an increasing need for technological development for efficient eRNA detection because of its physicochemical instability. The present study conducted a series of aquarium experiments using zebrafish (Danio rerio) and validated the methodologies for capture, preservation, and extraction of eRNA in a water sample. In the eRNA extraction experiment, an approximately 1.5-fold increase in lysis buffer volume resulted in a more than sixfold increase in target eRNA concentration. In the eRNA capture experiment, although GF/F and GF/A filters yielded similar eRNA concentrations, a GF/A filter may be capable of passing through more volume of water samples and consequently collecting more eRNA particles, given the time required for water filtration. In the eRNA preservation experiment, the use of RNA stabilization reagent (RNAlater) allowed for stably preserving target eRNA on a filter sample at - 20 and even 4 °C for 6 days at least. Altogether, the findings enable the improvement of eRNA availability from the field and easily preserve eRNA samples without deep-freezing, which will contribute to the refinement of eRNA analysis for biological and physiological monitoring in aquatic ecosystems.

Utilizing the state of environmental DNA (eDNA) to incorporate time-scale information into eDNA analysis.

Environmental DNA (eDNA) analysis allows cost-effective and non-destructive biomonitoring with a high detection sensitivity in terrestrial and aquatic environments. However, the eDNA results can sometimes include false-positive inferences of target organisms owing to the detection of aged eDNA that has long since been released from the individual and is more likely to be detected at a site further away from its source. In order to address the issue, this manuscript focuses on the state of eDNA, proposing new methodologies to estimate the age of eDNA: (1) DNA damage rate, (2) eDNA particle size distribution, and (3) viable cell-derived eDNA. In addition, the manuscript also focuses on the shorter persistence of environmental RNA (eRNA) compared with eDNA, highlighting the application of eRNA and environmental nucleic acid ratio for assessing the age of the genetic materials in water. Although substantial further research is essential to support the feasibility of these methodologies, incorporating time-scale information into eDNA analysis would update current eDNA analysis, improve the accuracy and reliability of eDNA-based monitoring, and further refine eDNA analysis as a useful monitoring tool in ecology, fisheries and various environmental sciences.

Can nuclear aquatic environmental DNA be a genetic marker for the accurate estimation of species abundance?

Environmental DNA (eDNA) analysis is a promising tool for the sensitive and effective monitoring of species distribution and abundance. Traditional eDNA analysis has targeted mitochondrial DNA (mtDNA) fragments due to their abundance in cells; however, the quantification may vary depending on cell type and physiology. Conversely, some recent eDNA studies have targeted multi-copy nuclear DNA (nuDNA) fragments, such as ribosomal RNA genes, in water, and reported a higher detectability and more rapid degradation than mitochondrial eDNA (mt-eDNA). These properties suggest that nuclear eDNA (nu-eDNA) may be useful for the accurate estimation of species abundance relative to mt-eDNA, but which remains unclear. In this study, we compiled previous studies and re-analyzed the relationships between mt- and nu-eDNA concentration and species abundance by comparing the R2 values of the linear regression. We then performed an aquarium experiment using zebrafish (Danio rerio) to compare the relationships across genetic regions, including single-copy nuDNA. We found more accurate relationships between multi-copy nu-eDNA and species abundance than mt-eDNA in these datasets, although the difference was not significant upon weighted-averaging the R2 values. Moreover, we compared the decay rate constants of zebrafish eDNA across genetic regions and found that multi-copy nu-eDNA degraded faster than mt-eDNA under pH 7, implying a quick turnover of multi-copy nu-eDNA in the field. Although further empirical studies of nu-eDNA applications are necessary to support our findings, this study provides the groundwork for improving the estimation accuracy of species abundance via eDNA analysis.

Complex interactions between environmental DNA (eDNA) state and water chemistries on eDNA persistence suggested by meta-analyses.

Understanding the processes of environmental DNA (eDNA) persistence and degradation is essential to determine the spatiotemporal scale of eDNA signals and accurately estimate species distribution. The effects of environmental factors on eDNA persistence have previously been examined; however, the influence of the physiochemical and molecular states of eDNA on its persistence is not completely understood. Here, we performed meta-analyses including 26 previously published papers on the estimation of first-order eDNA decay rate constants, and assessed the effects of filter pore size, DNA fragment size, target gene, and environmental conditions on eDNA decay rates. Almost all supported models included the interactions between the filter pore size and water temperature, between the target gene and water temperature, and between the target gene and water source, implying the influence of complex interactions between the eDNA state and environmental conditions on eDNA persistence. These findings were generally consistent with the results of a reanalysis of a previous tank experiment which measured the time-series changes in marine fish eDNA concentrations in multiple size fractions after fish removal. Our results suggest that the mechanism of eDNA persistence and degradation cannot be fully understood without knowing not only environmental factors but also cellular and molecular states of eDNA in water. Further verification of the relationship between eDNA state and persistence is required by obtaining more information on eDNA persistence in various experimental and environmental conditions, which will enhance our knowledge on eDNA persistence and support our findings.

Selective collection of long fragments of environmental DNA using larger pore size filter.

Environmental DNA (eDNA) can exist in water with various sizes and states. Among them, relative to extra-cellular DNA, intra-cellular DNA such as cell and tissue fragments can mainly be detected at larger size fractions, and may be protected from enzymatic DNA degradation processes. Here, we verified the hypothesis that the selective collection of such large-sized eDNA enhances the efficiency of capturing less-degraded eDNA, based on a tank experiment using Japanese Jack Mackerel (Trachurus japonicus) as a model species. We concentrated different volumes of rearing water using the filters with different pore sizes (0.7 µm and 2.7 µm), and quantified the copy number of short and long mitochondrial and short nuclear DNA fragments of target species in water samples. As a result, the ratio of long to short eDNA concentrations was higher in the larger pore size filter, which would support our stated hypothesis. In addition, the ratio of nuclear to mitochondrial eDNA was lower in the larger pore size filter. These results imply a difference in the persistence of nuclear and mitochondrial DNA between intra- and extra-cellular environments. Moreover, larger filter pore size did not necessarily decrease the yields of eDNA, and there was little difference in yields in smaller filtration volumes. The findings of this study indicate the potential to select information from eDNA signals by focusing on eDNA of specific size and state, which may contribute to efficient utilization of the information on species taxonomy and physiology in water samples.

Particle Size Distribution of Environmental DNA from the Nuclei of Marine Fish.

Environmental DNA (eDNA) analyses have enabled a more efficient surveillance of species distribution and composition than conventional methods. However, the characteristics and dynamics of eDNA (e.g., origin, state, transport, and fate) remain unknown. This is especially limited for the eDNA derived from nuclei (nu-eDNA), which has recently been used in eDNA analyses. Here, we compared the particle size distribution (PSD) of nu-eDNA from Japanese Jack Mackerel (Trachurus japonicus) with that of mt-eDNA (eDNA derived from mitochondria) reported in previous studies. We repeatedly sampled rearing water from the tanks under multiple temperatures and fish biomass levels, and quantified the copy numbers of size-fractioned nu-eDNA. We found that the concentration of nu-eDNA was higher than that of mt-eDNA at 3-10 µm size fraction. Moreover, at the 0.8-3 µm and 0.4-0.8 µm size fractions, eDNA concentrations of both types increased with higher temperature and their degradation tended to be suppressed. These results imply that the production of eDNA from large to small size fractions could buffer the degradation of small-sized eDNA, which could improve its persistence in water. Our findings will contribute to refine the difference between nu- and mt-eDNA properties, and assist eDNA analyses as an efficient tool for the conservation of aquatic species.

Effect of water temperature and fish biomass on environmental DNA shedding, degradation, and size distribution.

Environmental DNA (eDNA) analysis has successfully detected organisms in various aquatic environments. However, there is little basic information on eDNA, including the eDNA shedding and degradation processes. This study focused on water temperature and fish biomass and showed that eDNA shedding, degradation, and size distribution varied depending on water temperature and fish biomass. The tank experiments consisted of four temperature levels and three fish biomass levels. The total eDNA and size-fractioned eDNA from Japanese Jack Mackerels (Trachurus japonicus) were quantified before and after removing the fish. The results showed that the eDNA shedding rate increased at higher water temperature and larger fish biomass, and the eDNA decay rate also increased at higher temperature and fish biomass. In addition, the small-sized eDNA fractions were proportionally larger at higher temperatures, and these proportions varied among fish biomass. After removing the fish from the tanks, the percentage of eDNA temporally decreased when the eDNA size fraction was >10 µm, while the smaller size fractions increased. These results have the potential to make the use of eDNA analysis more widespread in the future.

Rapid degradation of longer DNA fragments enables the improved estimation of distribution and biomass using environmental DNA.

The advent of environmental DNA (eDNA) analysis methods has enabled rapid and wide-range ecological monitoring in aquatic ecosystems, but there is a dearth of information on eDNA degradation. The results of previous studies suggest that the decay rate of eDNA varies depending on the length of DNA fragments. To examine this hypothesis, we compared temporal change in copy number of long eDNA fragments (719 bp) with that of short eDNA fragments (127 bp). First, we isolated rearing water from a target fish species, Japanese Jack Mackerel (Trachurus japonicus), and then quantified the copy number of the long and short eDNA fragments in 1 L water samples after isolating the water from the fish. Long DNA fragments showed a higher decay rate than short fragments. Next, we measured the eDNA copy numbers of long and short DNA fragments using field samples, and compared them with fish biomass as measured by echo intensity. Although a previous study suggested that short eDNA fragments could be overestimated because of nontarget eDNA from a nearby fish market and carcasses, the eDNA concentrations of long fragments were correlated with echo intensity. This suggests that the concentration of longer eDNA fragments reflects fish biomass more accurately than the previous study by removing the effects of the fish market and carcasses. The length-related differences in eDNA have a substantial potential to improve estimation of species biomass.